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Abstract 

We provide a representation of the maximal difference between a stan- 
dard Brownian bridge and its concave majorant on the unit interval, from 
which we deduce expressions for the distribution and density functions 
and moments of this difference. This maximal difference has an applica- 
tion in nonparametric statistics where it arises in testing monotonicity of 
a density or regression curve. 



1 Introduction 

Motivated by applications to the theory of nonparametric statistics, indicated 
later in this introduction, we provide a useful representation of the maximal 
difference 

M := sup (\<3 [Qtl] B](u)-B(u)) 
ue[o,i] 

where (B(u),0 < u < 1) is a Brownian motion, and Cjf denotes the (least) 
concave majorant of a function / denned on an interval /. See Figure [TJ 

Our representation of M is presented in the following theorem, in terms of 
the distribution of M3, the maximum of a standard Brownian excursion, which 
can be represented also as the maximum of a three dimensional Besscl bridge, 
or as 

M 3 = sup B hl \u) - inf B hl \u) (1) 
«e[o,i] u6[o,i] 

where B hr denotes the standard Brownian bridge obained by conditioning B 
on B(l) = 0. See Q2], [5] for background and further information about the 
distribution of M3, which has been extensively studied. 
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Theorem 1 The distribution of the maximal difference M between a Brownian 
motion B and its concave majorant is determined by the identity in distribution 

M = max JL~M Z i (2) 

3 

where 

• (Lj,j = 1,2. ..) is the uniform stick- breaking process 

L 1 :=W ll L 2 := (1 - W 1 )W 2 , L 3 := (1 - Wi)[l - W 2 )W 3 , . . . 

derived from a sequence Wi,W 2 , . . ■ of independent uniform (0,1) vari- 
ables; 

• [M 3 j,j = 1,2...) is a sequence of independent random variables, each 
with the distribution of M%; 

• the two sequences (Wj) and (M 3 j) are independent. 

Moreover, M is independent of B(l), so the identity in distribution (|2|) holds also 
if the Brownian motion B on [0, 1] is replaced by the Brownian bridge obtained 
by conditioning B on B{1) = b, for arbitrary real b. 

To be more explicit, following Groeneboom [10], we observe that the vertices 
of the concave majorant ([C[ 0j i]B](w), < u < 1) partition [0, 1] into a countable 
collection of subintervals, with accumulation of vertices at both and 1, but 
with only a finite number of vertices almost surely in (e, 1 — e) for each e > 0. 
The sequence (Lj) represents the lengths of these maximal subintervals, over 
each of which the concave majorant has a segment with particular slope. These 
Lj arc arranged in a suitable random order, while y^LjM^j represents the 
maximum value of [C[ i^B](u) — B(u) for u in the interval of length Lj. Here 
and throughout the paper, the length of a segment of a concave majorant refers 
to the length of the time interval associated with the segment, rather than the 
length of the segment in a two-dimensional picture such as Figure [TJ 

For the proof of Theorem Q] we combine two different ingredients: 

• Groeneboom's description [10] of the concave majorant of Brownian mo- 
tion on the infinite interval [0,oo); 

• Suidan's description [25| of the joint law of ranked lengths of intervals in 
the partition of the time interval [0, 1] generated by vertices of C[ 0j i]B. 

In principle, the second of these ingredients must be derivable from the first. Fol- 
lowing Groeneboom, we use Doob's transformation to map Brownian motion on 
[0, oo) to Brownian bridge on [0, 1], and this mapping determines the law of the 
interval partition of [0, 1] derived from C[o,i]£? br . While wc provide some details 
of this in Section [4j we are unable to fully derive the stick-breaking representa- 
tion of interval lengths this way. Still, by developing the results of Groeneboom 
and Suidan, and by exploiting the fact that a uniform stick-breaking process is 
invariant under a size-biased random permutation (see [11], [18], and [20]) we 
see that the Lj appearing in Theorem [T] can be constructed as follows: 
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Figure 1: Plot of a Brownian bridge and its concave majorant. The bullets 
depict the vertices, whose times are indicated by the tick marks on the x-axis. 
The length of the dotted vertical segment gives the maximum difference between 
the Brownian bridge and the concave majorant. The Haar approximation was 
used to generate the Brownian bridge on a discrete partition of [0, 1] with a 
mesh equal to 2~ 15 . 

Corollary 2 Let Ui, 11%, . . . be a sequence of independent uniform (0, 1) vari- 
ables, independent of the Brownian motion B, and let L 1 ,L 2l --- be the size- 
biased random permutation of lengths of segments o/C[o.i]-B defined by: 

• L\ is the length of the interval containing U\, 

• L2 is the length of the interval containing the first point Ui that is not in 
this interval, 

• L3 is the length of the interval containing the next point Uj that is not in 
either of the first two intervals, 
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and so on. Then 

L^L^il-LtlLs/il-L, -L a ),... 

is a sequence of independent uniform (0, 1) random variables, and this sequence 
is independent of B{1). 

In particular. 

• the length L\ of the segment of C[ 0j i]i? br covering a uniform random 
number U\ is itself uniformly distributed in [0, 1]. 

We show in an appendix how this fact can be verified from Groencboom's [lOj 
joint density for the location of vertices of the concave majorant of Brownian 
motion on [0, oo). But the computations are difficult and we are currently 
unable to extend this method to obtain the assertions about Lj for j > 2 from 
Groeneboom's results. Rather, we rely completely on Suidan's approach for this 
part of the argument, which is essential for our proof of Theorem [1] 

To motivate our study of the distribution of M, we recall that in testing 
whether a density or a regression function on [0, 1] is decreasing, the supre- 
mum distance between the empirical estimator of the function and its concave 
majorant is used as a test statistic; see Kulikov and Lopuhaa, [TJ and Durot 
[5]. From |15| it follows that the the supremum distance between the empiri- 
cal distribution function and its concave majorant attains its maximum at the 
uniform density on [0,1]. Moreover, it is known that this statistic, when mul- 
tiplied by y/n, converges in distribution to the maximal difference M between 
a Brownian motion B on [0, 1] and its concave majorant. For the regression 
problem, the distribution of M appears again as the limit of a similar scaling of 
the supremum distance between the cumulative sum diagram and its concave 
majorant if the true decreasing regression is constant, which is also known to be 
the least favorable regression function for the testing problem in question [8]. 
Durot [8] established continuity of the distribution of M. It is an immediate 
corollary of Theorem [1] that M in fact has a density, and we provide formulas 
for this density and for moments of M in Section [3] Additionally we give an 
alternative characterization of M based on the inverse of the Laplace transform 
of a function involving modified Bcssel functions of the second kind. This other 
characterization suggests a way of calculating the quantiles of M at any desired 
precision using some appropriate approximation method for the inverse of a 
Laplace transform, but we do not pursue this here. 

2 The concave majorant of Brownian motion on 
a finite interval 

We need to show that conditionally given C[ 0j i]-B the difference process C[ 0i i]-B— 
B behaves like a succession of independent Brownian excursions between the 
vertex times of Cm^B. Grocncboom [TU] proved a corresponding result for a 
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Brownian motion on [0,oo), stated as Theorem [3] below. From this result on 
[0, oo), we will prove a similar theorem on [0, 1] by using a space-time transfor- 
mation. 

For a fixed a > let 

Z a := max{_B(t) — at} 

and D a the time at which the maximum is attained; see also [7|. The point 
(D a , Z a + aD a ) is one vertex of C[o j00 )B, and it follows from [TOl Theorem 2.1] 
that (D\/ a ,a > 0) is a pure jump process. Let So > Si > S% > ■ ■ ■ denote 
the successive slopes of Cro l00 )-B to the right of D a , and S-i < < ■■■ the 
successive slopes of C[ 0jOO )-B to the left of D a , so S < a < 5_i almost surely. 
For i £ Z let Tj denote the length of the interval on which the slope of the 
concave majorant is Si, so that {Vi := J2j<i i £ Z} is the sequence of times 
of vertices of C[ 0iO o)-B. 

Theorem 3 (Groeneboom \10\/ ) Conditionally given C[ 0oo )-B, the difference 
process (Cr 0jOO )i?(t) — B(t),t > 0), is a succession of independent Brownian 
excursions of prescribed lengths T{ = Vi — Vi-i, i £ Z between zeros at the times 
Vi,i £ Z. 

See [10], [17], [4] for various proofs of this result. 

It is important now to distinguish clearly the restriction to [0, 1] of C[o iOQ ^B, 
the concave majorant of a Brownian motion B on [0, oo), and CVii-B, the 
concave majorant of a Brownian motion B on [0, 1]. The former concave ma- 
jorant has vertices accumulating only at 0, whereas the latter, which is the 
subject of the following Proposition |4] has vertices accumulating at both and 
1. These two concave majorant processes agree on some random interval [0,i?] 
with < R < 1 almost surely, and then differ on (R, 1], where R is the time of 
the last vertex of C[ 0jOO )-B before time 1, whose distribution is determined by 
formula (|16p . See also the remark of Groeneboom [10] p. 1022]. 

Proposition 4 Let {Vi,i £ Z} denote an indexing of the times of vertices of 
C[ ,i]-B so that Vi is an increasing function of i £ Z with limj—j-oo Vi = 
and lim^oo Vi = 1. Moroever, we assume that this indexing depends only on 

C[Q,1]B. 

(i) The random set of times of vertices {Vi,i £ Z} is independent of B{1). 

(ii) The difference process (Cp i]B(u) — £?(«), < u < 1) is independent of 
5(1). ' 

(Hi) Conseguently, the distribution of both the random set of vertex times and 
of the difference process is the same for an unconditioned Brownian motion 
B as it is for a Brownian bridge from (0,0) to (1,6) for every real value 
ofb. 
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Proof. The first two assertions follow easily from the facts that 



• the set of times of vertices of the concave majorant of a function (/(it), < 
u < 1) is the same for /(it) as it is for /(it) + cu for arbitrary real c; 

• for a Brownian motion B, the process (B(u) — uB(l),0 < u < 1) is a 
standard Brownian bridge independent of B{\). 

The third assertion follows easily from the first two. □ 

Theorem 5 Let X be either an unconditioned Brownian motion, or a Brownian 
bridge from (0,0) to (1,6) for some real 6, and {Vi,i £ Z} an indexing of 
the times of vertices ofCmuX as in Proposition^ Then conditionally given 
Vi = Vi,i £ Z the difference process (Cr 01 iX(ii) — X(u),0 < u < 1) is a 
concatenation of independent Brownian excursions of prescribed lengths Vi—Vi—i 
between zeros at the times v^. 

Proof. Based on Proposition 0] (iii), it suffices to consider the case with 
B replaced by a standard Brownian bridge B hr . Then, according to Doob's 
transformation, 



It follows that the times Vi of vertices of the concave majorant of B hr on [0, 1] 
may be constructed as Vi = Ti/{Ti + 1) where the Tj's are the times of vertices of 
the concave majorant of B on [0, oo). Moreover, the difference between B hl and 
its concave majorant on (V^, Vi+i) is a transformation of the difference between 
B and its concave majorant on (Tj,Tj_|_x)- 

We will show now that this transformation maps Brownian excursions to 
Brownian excursions. Indeed, observe that the transformation between (it, £? br (it)) 
and (t,B(t)) is the restriction to a Brownian path of the space-time transfor- 
mation 



where 0<it<l, 0<t< oo, and x and y both range over all real numbers. 

Consider now the conditioning of B hl on B hY (u) = x and B hT {u) = x for 
some < u < u < 1 and real numbers x and x. The process 




where B is the standard Brownian motion 
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is a standard Brownian bridge in terms of which the path of B hr on [u, u] is 
represented as 

B hl '(u') = x + ^~^i,X* ( ^—^] +^—^(x-x), u'G [«,«]. 

\u — u J u — u 

On the other hand, with (i, y) = T(«, x) and (t, y) = T(u, i) the process 

r (w) := -j= , < u> < 1 

\Ji-t 

is another standard Brownian bridge in terms of which the path of B on [t, t\ is 
represented as 

B(t>) = y + Vi~t Y* (|f£) + t-±(y y), t> G [t, t\. 

It follows that for an arbitrary choice ofO<u<u<l and real values of x 
and x, with (t, y) and (t, y) the images of (u, x) and (u, x) respectively via the 
space-time transformation T in (|3|) , 

the image via the space-time transformation of a Brownian bridge 
from (u,x) to (u,x) is a Brownian bridge from (t,y) to (t,y). 

The key observation is that similarly, 

the image via the space-time transformation of the straight line from 
(w, x) to (u, x) minus a Brownian excursion of length u — u shifted 
to start at time u, is a straight line from (t, y) to (t, y) minus a 
Brownian excursion of length t — t, shifted to start at time t. 

Intuitively this is clear from the bridge result, by conditioning each of the bridges 
to stay above the line joining its endpoints. This can be made rigorous by a 
weak convergence argument, conditioning one of the bridges to go no more than 
a small distance e above the line, passing to the limit as e — > 0, and appealing 
to the result of [6]. □ 

The limiting argument at the end of the previous proof can also be reduced 
to an invariance of laws of standard Brownian bridges and excursions under a 
family of space-time transformations indexed by a pair of parameters < u < 
u < 1. To see this, observe that the relation between B hl and B implies that 
the standard bridge X* derived from B hr on [u, u] and the standard bridge Y* 
derived from B on [t,t\ are related by 

x» = l -"- t,(fl - u) y*f 1 o<v<i. (4) 

-u) V 1 -u-v(u-u)J' 

For the sake of completeness, we give the following lemma, which shows that the 
mapping from one standard bridge or excursion to another standard bridge or 
excursion which is implicit in Doob's space-time transformation is non-trivial, 
but nonetheless easily checked: 



7 



Lemma 6 For each fixed choice ofO<u<u<lif (Y*(v), < v < 1) is a 
standard Brownian bridge, then (X*(v),0 < v < 1) defined by has the same 
distribution as (Y*(v),0 < V < 1). Moreover, the same is true with the standard 
Brownian bridge replaced by a standard Brownian excursion, or by a standard 
Bessel bridge of any dimension. 

Proof. For the standard Brownian bridge, the result can be derived as above, 
or checked more directly by observing that Y* is evidently a centered Gaussian 
process with continuous paths, so it suffices to check that K(X* (v)X* (w)) = 
K(Y* (v)Y* (w)) = v(l—w) for < v < w < 1, and this is easily done. The result 
for a Bessel bridge of dimension d = 1,2,... follows easily from the representa- 
tion of the square of this process as the sum of squares of d independent standard 
Brownian bridges. For d = 3 this gives the result for standard excursion, by the 
well known identification of the standard excursion with a three-dimensional 
Bessel bridge due to David Williams. □ 

3 The distribution of M 

Proof of Theorem[TJ Proposition|4]implies that the distribution of the max- 
imal difference M is the same for a standard Brownian motion on [0, 1] as for 
a Brownian bridge from (0,0) to (1,6) for an arbitrary real number b. So it is 
enough to establish the characterization of M provided in Theorem [1] with a 
standard Brownian motion B on [0, 1] replaced by a standard Brownian bridge 
B hr . 

Let (Vi,i £ Z) denote the sequence of times of vertices of C[o i i]B br indexed 
in the same way as in Proposition [4] and Theorem and let (Tj, i £ Z) be the 
corresponding lengths of the segments of C[ ^]B hr ; i.e, Tj = V* — Vi-i,i £ Z so 
the sequence (T i: i £ Z) is independent of the sequence (M 3ji ,i £ Z). 

From the definition of M and Theorem [5j we have readily that 

M = max y/TiM^i 

i 

where the M$^s are independent random variables identically distributed as 
M3, the maximum of a standard Brownian excursion of length 1. 

Now let > T^ 2 ) > • • • denote the values of T, put in decreasing or- 
der. According to the result of Suidan [25j . the distribution of the sequence 
(TV',j = 1,2 . . .) is the limiting distribution of ranked lengths of cycles of a 
uniform random permutation of n elements, with cycle lengths normalized by n, 
commonly known as the Poisson-Dirichlet distribution with parameters and 
1. See also [9], [23], [19]. Suidan gives this result for the concave majorant of 
Brownian motion on [0, 1], but it applies just as well to the Brownian bridge, by 
Proposition 4. It is also known [26] that this asymptotic distribution of cycle 
lengths is obtained by ranking the terms of a uniform stick-breaking process in 
decreasing order. Theorem [T] and Corollary [2] follow immediately. □ 

Theorem [1] offers an easy way of simulating values of M. For the Monte 
Carlo implementation, we can use the representation |T]) of M3 and the Donskcr 
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Figure 2: In solid line is the plot of a Monte Carlo approximation of Fm based 
on a sample of size 5,000. In dashed line is the plot of the distribution function 
F3. As the figure suggests, F 3 (u) < Fm{u) for all u. This is an immediate 
consequence of ([T]) and the definition of M. 

approximation of B hl ' . It follows also from Theorem [1] that the distribution and 
density functions of M, denoted hereafter by Fm and Jm respectively, are given 

by 



F M {x) = E 



x > 



and 



/m(*)=X> 



1 f 1 



IPs 



x > 



(5) 



(6) 



where F% and fs are the distribution and density functions of M3, known to be 
given by 



F 3 (y) = l-2^(4 ?1 y-i)e- 2 "V y > 



(7) 
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and 

oo 

h(y) = 8 £ n 2 y(4n 2 y 2 ~ Z)e- 2n ^\ y > 0, 

71=1 

see e.g. |13j . [5] for further information about this distribution which has been 
extensively studied. Monte Carlo approximations of the distribution function 
Fm and its density /m based on 5,000 independent copies of the uniform stick 
breaking process are shown in Figure [5] and Figure [3] along with the distribution 
and density functions of M3 . 




Figure 3: In solid line is the plot of a Monte Carlo approximation of /m based 
on a sample of size 5,000. In dashed line is the plot of the density function ^3. 



Another technique for getting at features of the distribution of M is to 
introduce a standard exponential variable 71 independent of M, and observe 

that 

v /t7M = max \JiiLiM 3 ,i (8) 

i 

where the 71 Li are the points of a Poisson point process with intensity x~ 1 e~ x dx 
for x > 0, independent of the M 3ji 's. These are the jumps of a gamma process 
(7s j < s < 1) from which the Li can be recovered by first normalizing the 
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jumps by 71 and then putting the jumps in size-biased random order. This 
leads to the following proposition. 

Proposition 7 For all x > 

P( x /7lM <x) = exp(-v(x, 00)) 

where 



v\x, 00 



1-F* 



VvJ 



dy. (9) 

Proof. Note that the event {^/jiM < x} occurs if and only if there is no 
point \/"fiLiM^^ in the interval (x,oo). Since these y/jiLiM^^ are the points 
of a Poisson point process with intensity measure v; the image via the map 
(x, m) 1— > \fxm of the measure 

x- x e- x V(Ah e dm), 

it follows that 

P( % /7lAf < x) = exp(-v(x, 00)) 

with ^(a;,oo) given in ([9jl. 

See also [H [31 HT] for similar calculations. 



□ 



The following theorem gives another representation of Fm and /m- Let K m 
be the modified Bcsscl function of the second kind and order m € [0,oo), and 
consider the function 

00 

G(f) = Y[ exp |-4 2V2t nK^V^t n) - K (2V2t n) X , t > 0. (10) 

Also, let £ _1 be the operator of inverse Laplace transform. 
Theorem 8 

(i) For any real r > 0, we have 

2 



E(M r ) 



r(r/2) 



t r - L (l-G(t))dt. 



(11) 



(ii) For all x > the distribution and density functions Fm and Jm are given 
respectively by 



f m {x) = \c- x {G{Vt)/t) 



(12) 



f M (x) = ^\c- 1 (l-G(Vt)) 



(13) 
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Proof. It follows from (J7J and © that 

v(x,+oo) = 2J2 y-^^a^y- 1 - l)e~ 2n2x2y 1 dy 
n=l Jo 

= 2V t- 1 e - 2 " 2 - 2t (2^ 1 - 1)6-*" dt 
n=l Jo 

OO 

= 2j2(2A n (x) - B n {x)) 



where 



and 



so that 



A n (x)= £(e"* 1 /t 2 ) (2n 2 x 2 ) = 2V2nxK 1 (2V2nx) 
B n (x) = Icie^ 1 /t)l (2n 2 x 2 ) = 2K (2V2nx), 



(14) 
(15) 



OO 

exp(-i/(x, 00)) = exp|-4 2^2x nK x {2^2x n) - K (2V2x n) j = G(#). 

n=l 

The derivation of the expressions in (fT4")l and (fT5")) is deferred to the appendix. 
Now for r > 



£( 7l r/2 M 



(^l--F 7 r/2 Mr (x)) dx 

^1 — cxp(i/(a; 1 / r , 00))^ (ix 

(l - G(x 1/r )j dx = ( rx r - 1 (l-G(x))da 



The claim in (i) follows now from independence of the random variables 71 and 
M. 

To show (ii), we use again the independence of 71 and M . We can write 
nVTiM < x) = J F M e~ s ds. 
Using the change of variable t = s/^/x, we get 

J F M e ~ tX dt = cxp(-^(Va?, oo))/x. 



Thus, for all i > 



1 



[ -== ] = [C- 1 (exp{-v(y/Z,oo))/x)] (t) 



12 



Tabic 1: Approximation of the moments of M of order r = 1, 8. 



r 


E(M r ) 


X 


1 


0.999399 


0.997366 


2 


1.060258 


1.056803 


3 


1.195155 


1.190869 


4 


1.431334 


1.427101 


5 


1.819154 


1.816777 


6 


2.448679 


2.452149 


7 


3.481508 


3.499897 


8 


5.212503 


5.266828 



and the expression of Fm follows. The expression /m can be obtained immedi- 
ately by using known properties of the operator C. □ 

The identity in (|11|) gives a way of calculating the rth moment of M via 
numerical integration. We used this approach to compute E(M r ) for r = 1, ..,8, 
and the values are reported in the first column of Table [TJ In the second column 
are the corresponding empirical estimators for the same moments using 20,000 
simulated independent copies of M . On the other hand, the inverse Laplace 
transforms in the expressions of Fm and Jm given in Q12[) and p3[) can be 
approximated with very high precision using for example the Gaver-Stehfest 
algorithm. We refer to Q] for a detailed description of this method as well as 
a nice summary of other classical algorithms used to approximate inverse of 
Laplace transforms. However, implementation of the Gaver-Stehfest algorithm 
requires using a multiple precision software or some adequate arbitrary precision 
library. This will be pursued elsewhere. 

4 Complements 

Groeneboom's description of C\o,<x,)B implies that for each fixed t > the joint 
density of V t ~ , the time of the last vertex of C[ 0jOO )-B before time t, and V t + , 
the time of first vertex of C[o iOC )-B after time t is given by the formula 

{0<v 1 <t<v 2 } 

where t/+ = yl y >o and X and Z are independent standard normal variables. We 
show in an appendix (Proposition [10]) that this joint density can be presented 
more explicitly as 

fv-,v+(?i,v 2 ) = <V2 \ Vi)2 (V^F~ arctan G/^F)) ^ 



fvr v+ ( u i ' "2) = 7 ToTv 



E 



X 



Z 



vi ^v 2 - Vl 
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Using the time transformation t 1— ► t/(t + 1), the joint density of X M and Y u , 
the last and first times vertices of the concave majorant of B hl occurring before 
and after u € (0, 1), is given by 



fx u ,Y v (x,y) 



1 



y-x 



Tr(y-x) 2 VV x(l-y) 

x l{0<K<u<y<l}- 



arctan 



y-x 
x(l - y) 



(17) 



Now let L\ be the length of the segment of the concave majorant of a stan- 
dard Brownian bridge B hl on [0, 1] covering a uniform random number U\ inde- 
pendent of B hl , and let us verify using this formula that the distribution of L\ 
is also uniform on [0, 1], as shown already by Corollary [2j Using the expression 
in (jTTJ) , we find after some algebra that L x has density 



1-/ 



1 



1 



x(l — I — x) 

1 



arctan 



I 



7Tl Jo I V^l - 



arctan 



x(l-l- x) J 

k w 

du, 



dx 



y/u(l ~ U) 



using the change of variable u = x/(l — I) and putting k = \Jlj (1 — I). 
Now, 



«r/2 

/ 2d9 = tt, by the change of variable u = sin 2 (0), 
Jo 



r 1 du 

Jo y/u(l - u) 
Using the same change of variable we can write 



arctan 



k 



y/u(l-u) 



du = 



tt/2 



arctan 



arctan 



2k 



sin(20) 
2k 



5"' 



TT 

2 



sin(t) 
— cos(t) arctan 
cos 2 (£) 



sm(26)d9 

sin(t)dt 
b 



sin(t) 



o sin 2 (t)+b 2 
" /2 cos 2 (t) 



dt, with b = 2k 



o sin 2 (£)+& 2 

tt/2 



dt 



dt 



b 2 + lj tan 2 (t)+6 2 /(6 2 + l) 



where the last equality follows from the identity 1 / cos 2 (t) = 1 + tan 2 (t) . Put 
c = b 2 1 [b 2 + 1). Using the change of variable tan(f) = x, we get 



nir/2 



dt 



l tan 2 (£) + c 



dx 



1 



(1 + x 2 ){x 2 + c) 



1 + X 2 



dx 



14 



1 7T 



Now, 



and hence 

r.1 



l-c2 \y/c 



4/ 



arctan 



21 + 1 + P 



du = — 
2 



4/ 



i - 



2V~l fl + 1 



- i 



irVl 
l + Vl 



which implies that /li(0 = 1 for < I < 1. Thus we have an independent 
verification of the uniform distribution of L\ asserted by Corollary [2j But it is 
far from clear how to derive the further results of Corollary [2] by this approach. 
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5 Appendix 

We collect in this appendix some further identities and computations of an 
analytic kind which arise in connection with the material of the paper. 

Proposition 9 Let X and Y have the standard bivariate normal distribution 
with correlation E(XY) = p. Then 



E{X+Y+) = — y/1 - pi + P P(X > 0, Y > 0) (18) 

where 

P(X > 0, Y > 0) = i- (| + arctan(p/v/l-p2)) (19) 
which for p < can be rewritten as 



P(X > 0,Y > 0) = -— arctan( A /l -p 2 /p). (20) 

Proof. Roscnbaum [22l formula (5)] gives an expression for E(XYl{ X >h,Y>k}) 
in terms of the probability P(X > h,Y > k) and the standard normal density 
and distribution functions. In the present case h = k = and Roscnbaum's 
formula simplifies to (fTg)) . Formula (fl^)) is well known, see e.g. If p < 

the passage from (fTT)]) to is made via the trigonometric identity arctan(a;) + 
arctan(l/x) = -n/2, Vi<0. □ 
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Proposition 10 Let a,b > 0, and Z and W independent standard normal vari- 
ables. Then 



E 



Z+ 



w z 

a b 



1 (b (b 

— — arctan — 

2o7r V a \a 



where y+ = yl y >o- 



Proof. Let t := b/a, and X := Z and Y = (tW — Z)/yT +_g. S ince (X,Y) 
is a standard bivariate normal with correlation p = —l/\Jt 2 + 1 < where 
t = — \J\ — p 2 1 p > then it follows from Proposition [5] that 

E[Z+(tW - Z) + ] = - JL (Vl - P 2 - parctan(v/l - p 2 /p) 

= — (t — arctan(i)) 

and the conclusion follows after replacing t by 6/a. □ 



Proof of the identities (|14|) and (|15p . We start by showing (fT5|) . Letting 
z = 2v / 2^a;, we need to show that 

i r°° e -l/t 2 , 
Ko(z) = H o e —e-^dt. 



We have that 



OO —lit POO 1 o 

6 e~*/*dt = e- / ^e-^dt 



t Jot 



and putting u = (zt — 2)/(2v / £), it follows that t — (u + \J u 2 + 2z) 2 / z 2 and 
dt = 2tdu/y / u 2 + 2z. Hence 

oo 1/t /-oo u 2 

e" z l 'Ut = 2e~ z / d M . 



* J-oo \/u 2 + 2z 



It follows that 

i r e-v* 



/ 4 cit = e -« / 6 cfa = K {z) 



2 Jo * J-oo Vu 2 + 2z 

see e.g. Hunter [12]. To show (fl4)) . it is enough to show that 

e -i/t 



/ 

JO 



f 2 



Using the fact linx^o zK\[z) = 1 (sec e.g. Mechel [IB]), and noting that the 
integral on the left hand side is equal to 1 for z = 1 it is enough to show that 

2 I —e-* 2 ^dt=(zK 1 (z)y. 
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The calculations above imply that 



We conclude by using the well-known identity (zK n (z)' = —zK n -i(z) where 
K m is the modified Bcsscl function of the second kind of order m. □ 
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